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Abstract 

We develop number theoretic tools that allow to perform computations relevant 
for the quantum mechanics over finite fields of arbitrary, odd size, with the same 
speedup that is enjoyed by the Fast Fourier Transform. 
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Recent developments in non-perturbative string theory,the discovery of D-branes 
and their dynamics, revealed a new sector of heavy solitonic objects through which dis- 
tances below string scales can be probed in the weak coupling regime (non-relativistic 
limit) 0. It seems that the most fundamental of these solitonic objects, the DO branes, 
have low energy effective lagrangian the 1-dimensional reduction of a lOd supersymmetric 
YM system,the so called SUSY YM quantum mechanics. Indeed a stack of N DO branes 
has SU{N) SUSY YM quantum mechanics as its low energy effective lagrangian 0, |^ 
and the target space collective coordinates of the DO branes become N hj N hermitian 
matrices functions of time[YM gauge potentials] and their SUSY partners. This implies 
that a non- commutative geometry setting is emerging for the description of the dynamics 
of DO branes H]. 

All the above has been lifted to the level of a candidate for the M-theory, (the the- 
ory which presumably unifies all known string theories), the famous by now M(atrix) 
Theory Curiously enough,the above picture resembles the (many years old) SU{N) 
truncation [0] of the excitations of the supermembrane theory in 11 dimensions Q in anal- 
ogy with the bosonic membrane SU{N) truncation[P, |10| . 

More close to the DO picture comes the work of [jlO|, where the discretized membrane 
and its noncommutative geometry,Finite Quantum Mechanics, was introduced, as a con- 
sistent truncation of the bosonic membrane and its dynamical symmetry that of the area 
preserving diffeomorphism group. In this discretized version of the membrane the elemen- 
tary excitations- degrees of freedom were assumed to be one particle states living on the 
membrane, like particles in discrete phase space Z^r x Zat. A physical analog system of 
these elementary excitations was proposed to be the Quantum Hall effect of one electron 
on a magnetic lattice of rational magnetic flux per plaquette.The Hilbert space of these 
elementary excitations is finite dimensional and the quantum mechanics of linear quantum 
maps was further developed in refs|Tl|, |12|, 0. The SU (N) matrices of the YM quantum 
mechanics can be thought as coherent states of such elementary excitations. The difference 
with the above mentioned model of elementary excitations of the membrane is that time 
is also discrete and the motion of these excitations is typically random and chaotic- a 
fact which at the quantum level is translated into extended, random wave functions for 
typical eigenstates. 

Although we are far from a realistic scenario for the role of these elementary excitations 
for the quantum dynamics of the SUSY SU{N) Quantum Mechanics we believe that 
further technical developments are necessary in order to acquire better understanding of 
the situation. 

On a more mathematical side Finite Quantum Mechanics [FQM] has been developed 
so far using representation theory of the modular group S'L (2, Zat), the linear canonical 
transormation group of the elementary excitations,for values of iV,prime or powers of 
primes. In this note we treat the case of general odd integers A^, using prime decomposition 
and the Chinese remainder theorem for the modular group and its represent at ions. The 
case of integers N = 2^ and general integers will be dealt with elsewhere. 

An immediate, practical, consequence of our work is the possibility to extend the Fast 
Fourier transform for any odd to the metaplectic representation of the modular group 
SL{2,Zn) 

We now recall the basic features of FQM. 



The torus phase space has been the simplest prototype for studying classical and 
quantum chaos [14, 15, 1^, l?]- Discrete elements of SL{2, M), i.e. elements of the modular 
group SL{2, Z), are studied on discretizations of the torus with rational coordinates of the 
same denominator /, {q,p) = {ni/l,n2/l) G F, ni,n2,l G Z and their periodic trajectories 
mod 1 are examined studying the periods of elements A G SL{2, Z) mod /. The action 
mod 1 becomes mod / on an equivalent torus, (ni, 712) G IT. The classical motion of such 



discrete dynamical systems is usually "maximally" disconnected and chaotic [T^, |T7|. 

FQM is the quantization of these discrete linear maps and the corresponding one-time- 
step evolution operators U{A) are / x / unitary matrices called quantum maps. In the 
literature |T^, these maps are determined semi-classically. In ref. |Tl| , |12[ the exact 
quantization of SL{2, Fp), where Fp is the simplest finite field of p elements with p a prime 
number was studied in detail. In ref. these results were extended to powers of primes, 
p"-the group is then SL{2, Zpu). 

The Hilbert space Hr of the wave functions on the torus F = C/L of complex modulus 
T = Ti + iT2, where L is the integer lattice, L = {mi + rm2|(mi, 7712) G Z x Z} is defined 
as the space of functions of complex argument z = x + iy 



iirn T+2-iTinz 



with norm 11 



e-'''y"/^'\fiz)\^dxdy, T2>0 (2) 
Consider the subspace Mi(T) of Hp with periodic Fourier coefficients {cn\n& of period / 

Cn = Cn+i n G Z, / G N (3) 
The space ]H[/(F) is /-dimensional and there is a discrete Heisenberg group |T^, with 



generators Si/i and Ti acting as |T9|, ^ 

iSl/lfKz) : 



27rin// 27rin.3+7rin^7 



(4) 



On the /-dimensional subspace of vectors (ci, . . . , q) the two generators are represented 
by 

('^l//)ni,n2 Qni,n2 ^^ni,n2 



('^)ni,n2 Pn\,n2 ^n\ — l,r, 

with uj = exp(27ri//). The Weyl relation becomes 



(5) 



QP = ooPQ 



(6) 



and the Heisenberg group elements are 

The generators Jr^s satisfy the following composition law 

Jr,sJr',s' = '^^^ '^'''^^^^ _^ ^'^ S + s' 

and the "commutation" relations 

The metaplectic representation of SL(2, Z;) is defined by the relation 

U-\A)Jr,sU{A) = Ji^r,s)A (10) 

where A is an element of SL{2,Zi). In the literature the metaplectic representation of 
SL{2, Zi), (the group of 2 x 2, integer valued matrices mod /), is known for / = ||21|| p| 



The Weyl-Fourier form of U{A) is [13 



where 



p"— 1 

U{A) = ^^^^ y ^^[br-H'i-ays-cs-mj^ ^ ^^^^ 
^ r,s=0 



p"-i 

= ^E^^^^ (12) 



0" X 



r=0 



All the operations in the exponent are carried out mod p". If 5 = mod we use 
the trick 

a6\_/0 \\(-c-d 
cd)~\-\^)\a b 



(13) 



and the fact that U{A) is a representation (cf. ref. [|1^] and below). 

We shall now work out some technical details of the representation theory of SL{2, I^n), 
that we need for the factorization of the Heisenberg group and of the metaplectic repre- 
sentation. 



We start with the Chinese remainder theorem for numbers ||23||. Let G Z be a 
non-prime, that may be written as a product of two co-prime factors, Ni and A^2, viz. 
N = N1N2. We denote by A^Z the set of all multiples of A^. Then any r G Z/A^Z may be 
written uniquely as 

r = rimirii + r2m2n2 



^The representation theory of the symplectic group SL{2, Fpn), over the finite field Fpn . may be found 
in ref. PI. 



where ri = rmodA'"i, r2 = rmodA^2, mi = N/Ni,m2 = N/N2, rii = m^^modNi,n2 — 
m2^modN2. 

In other words, we may estabhsh a one-to-one correspondance between the number 
r and the 2-tuple (ri,r2). This last defines the Sino representation of r. This 2-tuple 
may be promoted to a bona fide element of a set whose elements have the following 
properties, for any two numbers r, r' e Z/NZ: 

r X r' (rir^modA^i, r2r2modA^2) 

and 



r + r' (ri + r^modA^i, r2 -|- r2modA^2) 

It is immediate to generalize this result to the case where N = Ni x N2 x ■ ■ ■ where 
all pairs of factors are co-prime. The decomposition reads 

r = r-jriini -\- r2m2n2 H h rkUikTik 

where rrii = N/Ni,ni = m^^modNi and one may similarly establish a one-to-one corre- 
spondance between r and the /c— tuple (ri, r2, . . . , r^), element of the set Furthermore 
note that V^v has the property 

Using these relations it is now possible to establish that 

SL{2, Zn) = SL{2, ZtvJ x SL{2, Z^vJ x • • • x SL{2, ZnJ (14) 

Indeed, consider the case k — 2 and an element of 5'L(2, Zjv) of the form 

^^fa b \f (a„a2) (^biM) ^ 
\ c d J \ (ci,C2) (di,d2) J 

It will now be shown that this element is an element of the set SL[2,'Zn^) x SL{2,'Ln^) 
Consider a generic element of SL{2, ZatJ. It may be written as 

?A) (16) 



(ci,0) (di,l) 

Take now a generic element of 5'L(2, Zjvj), that may be written 

/ (l,a2) (0,62) \ (.r,. 

-^'-[{0,02) {l,d2)) (1^) 

It is straightforward to check that Ai ■ A2 = A. Using the Chinese remainder theorem 
for numbers we know that the decomposition is unique. 



Let us close with the remark that the matrices of the form 

a h 
—b a 



with + 6^ = ImodA^ generate the group 02{N) <l SL{2, Tjn)- Once more it is possible 
to show that 

02{N) = 02{Ni) X 02{N2) (19) 



VIZ. 

a h\^( (ai,l) (61, 0) \ ( (l,a2) (0,62) 
-ha I (-61,0) (ai,l) i I (0,-62) (1,^2) 



(20) 



In the following we discuss the implications of the factorization of SL{2, Z^v) (cf. previous 
discussion) for the metaplectic representation. We begin with the factorization of the 
Heisenberg group i){N = N1N2) (cf. the work of J. Schwinger in [[18| ). 

Indeed, using the Chinese remainder theorem and the commutation relations of the 
i7r,s, we find 

>Jr,s = J^nmini + r2m2^2, Simirii + S2m2n2 ~ 

JrimiTii, SiminiJr2'm2n2, S2Tn2n2 = 

>Jr2m2n2, S2?7i2^2^'"i^i^i, Simirii 
since the extra phase factor equals unity. 

We shall now use the factorization properties of the Heisenberg group generators, J7'r,s 
to obtain the decomposition of the unitary operator U{A), A G SL{2, Z^r), that governs 
the time evolution of the quantum system in the case in hand. To do this we recall that 
the evolution of the generators X,s is given by 

U-\A)Jr,sU{'^) = Jir,s)A (21) 

If iV = N1N2, then A = Ai- A2, with Ai G SL{2, ZjvJ and A2 G SL{2, Z^J. We shall 
show that U{A) = U{Ai) ■ t/(A). 

Proof: U {A) may be written as a linear combination of the generators Jr^s as 



U{A) = '-^^ ^ Jhr^ + {d-a)rs-cs-)/{28)j^^^ (22) 



r, s = 

where 

A = 



a h 
c d 

Using the Sino representation, xr"^ = (xi, a;2)(r^, r|) = (xirJ,X2r|) = Xir^nimi+X2r|n2m2 
and the double sum is seen to split into the product of two sums 



N-l Ni-l,N2-l 

r=0 ri,r2=0 l'^'X\ 

_L ^ g27rimiXirf/Ari j ^ _L ^ g27rim2X2r^/A/'2 



J_ Y^M^^^ = \^ p2T[\{xxr\nx'mx + X2rln2m2) / {N1N2 

N 4^ N1N2 



Ni ^ \N2 

^ ri=0 / \ ^ ri=0 



which leads to the relation 

a{x) = (7{miXi)cr{m2X2) (24) 
This takes care of the prefactor. The phase is re-arranged as follows: 

(}) - (&r2 + {d- a)rs - cs^)/{2S) = 
{{birl + {di - ai)riSi - CiSi)/(2(5i)) mini + 
((&2ri + {d2 - a2)r2S2 - C2S2) / {262)) m2n2 = 
(f)iN2mi + (j)2Nim2 



The upshot of this is that U{g) may be rewritten as 

AM 



J'l,?'2=0 

/ J ^ r\m\n\,s\m\n\ / ^ ^ N2 >-^r-2m2n2,S2m2n2 K'^^J 

^ ri=0 ^ r2=0 

As a consequence U {Ai) and (^2) commute. Now we establish an isomorphism between 
U{Ai) and U{A2) with Ui{Ai) and C/2(^2), where Ui and C/2 are the metaplectic repre- 
sentations of dimension Ni and N2 respectively. Indeed we shall exhibit a permutation 
matrix R with the properties 

RPR^ = Pi P2 
RQR^ ^Q^®Q2 



RU{g)R' =U,{g,)®U2{g2) 
where P and Q {Pi, Qi,i — 1,2) are the generators of the Heisenberg group '^{N) (resp. 

It is enough to prove that the matrix R has the property 

Rek = efci ® efe2 (25) 

where (e^.,-? = 1,2) are the eigenvectors of P (rcsp. Pj,i = 1,2) and the indices run 
as /c = 1, . . . , N,ki = 1, . . . , iVj,i = 1, 2. The other properties indeed are consequences of 
this. 

We start with an exphcit form for the eigenvectors of P 

^^k{l-l) 

ek,i^^=^, k,l^l,...,N (26) 
V iV 



Using the Sino representation we find 



ki{ji-l)mini+k2{j2-l)m2n2 , ki{ji-l)mi ^ k2{j2-l)m2 

(27) 



In order to construct the matrix R, we compare the rhs of eqs.(25,2^. The indices ji 
and j2, in eq.(^), from the Sino decomposition of j, run from 1 to A^^i (resp. A^2) and 
the corresponding values of j fill up a A^i x N2 array. On the other hand, the rhs of 
eq.(p5D defines, through the tensor product, a decomposition of j into indices ji and 
and defines another A^^i x N2 array, that is related to the previous one by a permutation 
matrix, namely R. Specifically, if we denote the Sino decomposition array by 

{Ji,j2} (28) 

the matrix R, with rows indexed by i and columns indexed by j, has elements equal to 1 
when {71,72} = and zero otherwise. This construction is straightforwardly generalized to 
more than two co-prime factors of (cf. the book by Schroeder in |Q). The stationary 



eigenvalue problem for the unitary evolution operator U is reduced to that corresponding 
to each individual factor of the tensor decomposition of the matrix U (e.g. in the case 
where Ni and N2 are powers of primes of the type 4/c + 1 explicit expressions are known 
for the eigenvectors and eigenvalues p^, |13|)- 

The matrix R is used in (classical) Fast Fourier algorithms to reduce the number of 
operations from 0{N^) to 0{N\ogN)^). 

By construction, therefore, the property in eq.(p5D holds and this implies the decom- 
position of the P operator. The decomposition of the Q operator follows immediately, 
since it is diagonal in this basis and the operations may be carried out element by element. 
From these both the decomposition of the J^r s and U{g) follow since (as may be checked) 
RR^ = I. 

We close with the following remarks. The case N = 2^ cannot be studied by the 
methods developed here and new ideas are required. This case is, of course, particularly 
interesting for computational reasons. Indeed, all existing Fast Fourier algorithms are 
given for N = 2^. For powers of primes a similar speedup of operations may also be 
obtained (cf. Schroeder in PSI). On the other hand, what we have achieved here is 



the construction of fast algorithms for any odd N and for any quantum map that is a 
metaplectic representation of SL{2,Zn). 

It would be interesting to implement such maps by quantum gates, as already proposed 
for the quantum Fourier transform . 
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